Deviations from Fermi-liquid behavior above Tc 
in 2D short coherence length superconductors 
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We show that there are quahtative differences between the temperature dependence of the spin and charge 
correlations in the normal state of the 2D attractive Hubbard model using quantum Monte Carlo simulations. 
The one-particle density of states shows a pseudogap above Tc with a depleted A'^(O) with decreasing T. The 
susceptibility Xa and the low frequency spin spectral weight track A''(0), which explains the spin-gap scaling: 
1/TiT ~ Xa{T)- However the charge channel is dominated by collective behavior and the compressibility 
dn/dn is T-independent. This anomalous "spin-charge separation" is shown to exist even at intermediate 
\U\ where the momentum distribution n(k) gives evidence for degenerate Fermi system. 

PACS numbers: 71.27.+a, 74.72.-h, 74.25.-q, 74.40. +k 



The unusual normal state properties of high Tc su- 
perconductors (SC) have led to many studies exploring 
possible scenarios for the breakdown of Fermi liquid the- 
ory (FLT) [|l|. Much of this effort has been in one of 
two directions: search for genuinely new non- Fermi liq- 
uid ground states arising from strong correlations in the 
proximity of an insulating state, or search for a low en- 
ergy scale such that for temperatures above it one would 
see deviations from FLT. In this paper we take a different 
point of view and ask the question: Is the normal state of 
a short coherence length SC necessarily a Fermi liquid? 
More specifically, if the ground state is a condensate of 
fermion pairs, and the phase transition leads to a degen- 
erate Fermi system above Tc, do the correlations in that 
normal state have to be described by Landau's FLT? If 
not, what are the characteristic deviations from FLT? 

For a weak coupling SC, with the coherence length 
^0 ^ a, the lattice spacing, the normal state is a FL. As 
the coherence length decreases, and Co/a is not a small 
parameter, the Tc at which phase coherence is established 
and the pair formation scale separate [pP]. Certainly in 
the opposite extreme of tightly bound pairs, the state 
above Tc is a normal Bose liquid. The question then is 
whether there is a broad intermediate coupling regime, 
especially in 2D, where the normal state has a "Fermi 
surface" and yet exhibits non-Fermi liquid correlations. 

The simplest lattice model within which this problem 
can be studied is the attractive Hubbard model where 
the pair-size can be tuned by varying the strength of 
the interaction. In the absence of a small parameter we 
use quantum Monte Carlo (MC) simulations [Q,|| to gain 
insight into the intermediate coupling regime above Tc- 
While the U < Q Hubbard model is not a realistic mi- 
croscopic model for the high Tc materials, it has two 
merits: its simplicity (fewest parameters) and the relia- 
bility of low temperature MC simulations, since it is one 
of the few interacting fermion model which is free of the 
"sign problem" at all densities M . 



We begin by summarizing our main results which are 
obtained in the normal state (T > Tc): 

(1) The momentum distribution n(k) shows structure 
reminiscent of a Fermi surface, though broadened by both 
thermal and interaction effects. 

(2) A pseudogap opens up in the one-particle den- 
sity of states (DOS) well above Tc. The DOS is 
strongly T-dependent - hence the notation Nt{0) - with 
dNT{0)/dT > 0. We obtain 7Vt(0) using a new method 
which avoids numerical analytic continuation. 

(3) At low T, the spin susceptibility Xs{T) = 7Vt(0), 
while the compressibility (dn/dfi) is T-independent. The 
spin and charge responses are thus qualitatively different. 

(4) The low frequency spectral weight K{q) = 
Imx(q, a;)/a; in the spin channel is of the form K{q) ~ 
Nt {0)/rQ{q) where the T-dependence is largely in the 
DOS and ro(q) is essentially the same as in non- 
interacting system. This naturally explains the spin-gap 
scaling 1/TiT = J2n^i^) '^ XsiT), noted in our pre- 
vious work ||], and observed in the underdoped high Tc 
cuprates [Q. 

These results show that the normal state of a short 
coherence length 2D SC exhibits marked deviations from 
usual Fermi liquid behavior. One obtains a kind of "spin- 
charge separation" with the spin correlations determined 
by one-particle excitations, while the charge channel is 
dominated by collective excitations. 

Consider the attractive Hubbard Hamiltonian H — 

-t lli,j;a cLcjT +h.c. -\U\Y,i ni]n^l+^lY.^■,a ^i-^^ ^hcrc 
the hopping is between nearest neighbor sites. All ener- 
gies are measured in units of t, and the lattice spacing 
a = 1. 

We first establish that in the regime |t/| = 4 and {n) = 
0.5 S that we focus on here, our finite system results 
are indeed characteristic of the normal state. For this 
we estimate the correlation length £,{T) from the spatial 
decay of the SC order parameter correlation function on 
systems of size up to L = 16. We find that for \U\ — 4 



we have ^-3-4atT= 1/6; thus ^(T) < L and the 
system is normal. (Details of this analysis and finite size 
scaling will be presented elsewhere). 

Single-particle properties: 

To ascertain that the system is degenerate at mod- 
erate \U\ we study the momentum distribution n(k) ~ 
^^(CkcrCko-) on large lattices of size up to iV = L^ = 16^; 
see Fig. 1. It is clear that n(k) shows a rapid variation 
with k - a "Fermi surface" - even though it is broad- 
ened by the temperature and also affected by the inter- 
actions. We emphasize that we are not in the large-|C/| 
preformed boson limit, where the constituent fermions 
(tightly bound into singlet pairs) are non-degenerate and 
their n(k) is independent of k. 

We now turn to the single-particle density of states 
(DOS) N{uj) with Lo measured from fi. This is given by 
N{uj) — j^ ^j^yl(k, w). Here A is the spectral function, 
which is related to the imaginary time Green function 
G(k,T) = -(T[ck.(T)cL(0)])via 



G(k,T)=- 



du ^"P^7"j , ^(k,.). (1) 



for < T < (3 — 1/T. To estimate A and from it the 
DOS, given MC data for G'(k, r), involves inverting this 
Laplace transform. We avoid the numerical complica- 
tions inherent in such an analytic continuation by deriv- 
ing a general expression for N(0) in terms of G'(k, r), 
which is valid provided there is no low energy scale in 
the problem. Fourier transforming to real space and 
looking at the local Green function at r = /3/2 we get 
G(r = 0,/3/2) = - J dLosech{Pu)/2)N{Lo)/2. Let n be 
the frequency scale on which there is structure in the 
DOS: for T < f7 we obtain [1 



iV(0) ~ -/3G(r = 0, r == /3/2)/7r. 



(2) 



The one-particle DOS obtained from (g) is plotted in 
Fig. 2 as a function of temperature. We see that for 
T > Tc a pseudo-gap develops at the chemical potential 
and the DOS is depleted as T is reduced. To empha- 
size this T-dependence of the DOS we use the notation 
Nt{0). This behavior should be compared with that in 
weak coupling where the DOS remains featureless above 
Tc, except for a small fluctuation dip [|0|. The pseudogap 
at intermediate coupling may be thought of as the evo- 
lution of the weak coupling fluctuation dip into a regime 
where a/S,o is no longer a small parameter. 

Spin and charge correlations: In Fig. 2 we also 
plot the uniform, static spin susceptibility Xs- The T- 
dependcncc of Xs was already noted in ref . Iq] ; what we 
see here is that this T-dependence comes entirely from 
that of the one-particle DOS, so that Xs{T) = Nt{0) (to 
within the errors inherent in extracting the latter). 

It is interesting to ask whether the charge channel 
also exhibits the same pseudo-gap. The compressibility 
{dn/dfi) was obtained by numerically differentiating |l|] 



the average density (n) determined as a function of fi. 
We found that (dn/d/ji) shows significant system size de- 
pendence (much more than, e.g., Xs); small system data 
show (finite size) upturns which are pushed down to lower 
T with increasing L. The results on the largest lattice 
(L = 16) plotted in Fig. 3 show that the system becomes 
more compressible with increasing \U\ (see below). More 
significantly, in sharp contrast to the one-particle DOS, 
{dnldfi) is very weakly T-dependent. 

Within a simple RPA (p-h bubbles) we obtain 
{dn/dfi)^P'^ = 2No{0)/ [1 - \U\NoiO)]. This is valid only 
\U\No{0) -C 1, where it correctly explains the trend that 
attractive interactions increase (dn/dfi); for large \U\ 
RPA fails in that it predicts an entirely spurious instabil- 
ity (phase separation) when \U\No{0) = 1. In fact, pairs 
form at large \U\ and their residual interactions are re- 
pulsive, so that the compressibility of the system remains 
finite for all \U\. It is simplest to see this in the broken 
symmetry state at T = where we find |12[| that, within 
mean field (MF) theory, {dn/dn) decreases monotonically 
with \U\ from {dn/dii) ~ 2No{0) [1 + \U\No{0)] for small- 
\U\ to {dn/dii) ~ \U\lAdt^ - 2/\U\ for \U\/t > 1. The 
T — Q MF result ||l^] is also shown in Fig. 3 and is found 
to be of the right order-of-magnitude as the normal state 
MC result; note that we do not expect {dn/dji) to change 
dramatically as T goes through Tc- 

The difference between the spin and charge response 
functions could be characterized as a sort of spin charge 
separation p3| . In its mildest form this exists even in 
a Landau Fermi liquid where the Xs a-nd (dn/dfi) are 
quantitatively different, the two being renormalized by 
different FL parameters: Fq and Fq. What we see 
here is much stronger: as a result of strong interac- 
tions, Xs and (dn/dfi) acquire qualitatively different T- 
dependences. As argued above, the spin response is dom- 
inated by incoherent single-particle excitations, which is 
T-dependent because triplet excitations require breaking 
up the singlet pair correlations while the charge channel 
is dominated by collective behavior. 

Spin-gap scaling: We next use our results for the T- 
depcndcnt DOS Nt{0) to gain insight into the suppres- 
sion of low frequency spectral weight in the spin channel 
as probed by -ftr(q) = hmjj^oInix(q, a')/a;. To contrast 
with our MC results it may be useful to recall that for 
a Fermi liquid (all quantities denoted by a subscript 0): 
A'o(q) = A^o(0)/ro(q) is T- independent for T < e^ with 
ro(q) r^ vpq- Further EqToHq) - ^o(O) leads to the 
Korringa law (l/TiT)o = Eq^o(q) - ^o(O). 

In Fig. 4 we plot the MC results for A'(q) for q ^ 
The analytic continuation required to obtain K was done 
using the method of ref. |^. From Fig. 4(a) we see that 
Ar(q) is more or less uniformly suppressed at all q with 
decreasing temperature. Further this T-dependence is 
similar to that of Xs{T) (or the DOS) shown at q = 
in the Fig. 4(a). For fixed T the q-dependence of K{q) 



resembles that of the non-interacting system as shown in 
in Fig. 4(b). 

The numerical results thus suggest that K{q;T) = ax 
A'T(0)/ro(q), with a independent of T and q, i.e., the T- 
dependence comes from a DOS with a pseudo-gap while 
the q-dependence is that of the non-interacting system. 
This form for K{q;T) leads to 1/TiT = Eq-f^(q;r) - 
aNo{0)NT{0), thus providing a natural explanation for 
the spin gap scaling 1/TiT ^ Xs{T) found in our earlier 
MC studies [§. 

We note that 1/TiT - Xs{T) is observed in the (pla- 
nar) O and Y NMR in a large number of (usually un- 
derdoped) cuprates 0. What distinguishes the results 
presented here from spin gap theories [|l4| based on spin 
models, is that the anomalies exist in a system with itin- 
erant carriers with a large Fermi surface. Secondly in 
the present approach these anomalies are directly related 
to the T-dependencc of the one-particle DOS. On the 
other hand, since we work with a "coarse-grained" model, 
rather than a realistic microscopic one, we cannot discuss 
the q-structure of the spin response leading to differences 
in the Cu and O NMR, or the doping dependences. It 
is also an important open problem to study the normal 
state precursors of short coherence length SC with gap 
anisotropy and/or nodes. 

In conclusion we have shown that the normal state of a 
2D short coherence length SC exhibits characteristic de- 
viations from Fermi liquid behavior: while the momen- 
tum distribution gives clear evidence for a degenerate 
Fermi system, the one-particle DOS shows a pseudo-gap 
p5[ and the spin and charge correlations show qualita- 
tively different temperature dependences. 
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FIG. 1. The momentum distribution n(k) versus [k| along 
[1, 0] and [1, 1] for \U\ = 4, (n) = 0.5 and T = 0.25. The Fermi 
function (U — 0) is shown as the long dashed curve and the 
U — —oo Bose limit result is plotted as the short dashed 
curve. The statistical error bars on all the MC data, unless 
explicitly shown, are less than the size of the symbols in this 
and other figures. 



FIG. 2. The one-particle density of states Nt{0) at the 
chemical potential (full triangles) and the spin susceptibility 
Xs (open squares) as a function of temperature for \U\ =4, 
(n) = 0.5 and L^ = 8^ 



FIG. 3. The compressibility (dn/d^i) for (n) = 0.5, |I7| = 
(full line) and |[/| = 4 (open circles), as a function of T ob- 
tained on a 16^ lattice. The T — non-interacting and the 
T — mean field result for |[/| = 4 are also shown. 



FIG. 4. (a) Low frequency spectral weight in the spin 
channel Jsr(q; T) — linitj^o Imx(q, uj)/uj for q 7^ (open sym- 
bols) plotted along [1,0] and [1, 1] for various T. The dashed 
lines are guides to the eye. The filled symbols plotted at q = 
show 2.0 X XaiT) with T corresponding to that of the open 
symbols. Note that the T-dependence of K{q) is similar to 
that of the susceptibility Xs{T). All of the data is for \U\ =4, 
(n) — 0.5 and L^ =8^. (b) At a fixed T the q-dependence 
of K{q) is qualitatively similar to that of the noninteracting 
case. To show this we compare the K{q) a.t T — 0.25 (open 
squares) with the full curve given by 2xs(r)7^°(q), where 
K^{q) is the essentially T- independent spectral weight for the 
non-interacting system. 
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